A method for characterizing coherent-state quantum gates 
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We discuss and implement experimentally a method for characterizing quantum gates operating 
on superpositions of coherent states. The peculiarity of this encoding of qubits is to work with 
a non-orthogonal basis, and therefore some technical limitations prevent us from using standard 
methods, such as process tomography We adopt a different technique, that relies on some a-priori 
knowledge about the physics underlying the functioning of the device. A parameter characterizing 
the global quality of the quantum gate is obtained by "virtually" processing an entangled state. 



Future quantum computers will be able to incorporate 
quantum logic in the treatment of information. For this 
purpose, operations are performed on quantum systems 
in order to implement the desired processing steps. Un- 
avoidably, such operations will only be an approximation 
of the ideal quantum gates, and the degree within which 
this approximation is acceptable is determined by the 
effectiveness of error correction codes [I] . The character- 
ization of quantum operations is then an important step 
to establish practical limits in the use of such devices. 

A crucial requirement is that the system remains under 
controlled manipulation, and is well-preserved from cou- 
pling to the environment, so that the action of the gates 
is not spoiled by decoherence. Optics offers an interesting 
option in this sense; within quantum optics, several pro- 
posals have been put forward, either based on a discrete 
[2], or a continuous variable approach [3]. 

Recently, a different point of view has emerged, which 
aims to join the strengths from both worlds, and en- 
codes quantum bits (qubits) in superpositions of weak 
coherent states , . The peculiarity of this approach 
is that the logical states |0) and |1) are represented by 
two non-orthogonal (thus non mutually-exclusive) states 
of the system, viz. two coherent states with the same 
amplitude and opposite phase \a) and \—a); this has 
demonstrated to be an error-correctable approach for 
moderately low intensities \a\ ^ 1.5, where the overlap is 
a)p ~ 10~^ [6 . The cost of this choice is that the 
gates are necessarily probabilistic, with a success rate 
that depends on this overlap; remarkably, though, the 
comparison with other optical schemes in terms of re- 
sources needed for scalability seems favorable [6, 7 . 

Gates in a coherent states architecture cannot be de- 
scribed by ordinary quantum operators, as they are de- 
fined using non-orthogonal states as a basis. More- 
over, decoherence processes most likely will not leave 
the output state inside the original reduced Hilbert space 
spanned by \a) and \—a). The current technologies do 
not even warrant that input test states belong to this 
space. Therefore, we cannot simply use standard tech- 
niques, such as process tomography [S^-HT, for a charac- 
terization. 

In this article we propose a way for characterizing these 
gates, which does not rely on a black-box approach [9^, 



but requires some modeling of the functioning of the gate. 
This is a realistic approach, since such a knowledge is 
needed to achieve the desired level of control. One can 
then identify a small number of parameters, accessible 
to the experimentalist, by which the gate process can be 
modeled. 

Our demonstration concerns the 7r-phase gate for 
coherent-state qubits proposed by Marek and Fiurasek 
[12] . The operating principle of the gate is that coherent 
states are eigenstates of the photon subtraction operator: 
a\ ±a) = ±a\ ± a). This operation thus corresponds to 
a TT-phase shift, up to an overall constant. In the lab- 
oratory, photon subtraction can be approximated by a 
beam splitter with low reflectivity, followed by detection 
on a photon counter. There have been demonstrations of 
such an effect [TBHlHl HQ. , but so far the characterization 
has focused on the states that could be produced by this 
technique, rather than on the device itself. While it has 
not been investigated if these gates are sufficient for full 
scalability, they nevertheless provide a reliable way to im- 
plement logical operations and test this architecture, in 
the same vein as other optical realizations [T3[|T4]. There- 
fore, they are an interesting testbed for our method. 

Testing quantum gates demands to operate them with 
at least two different orthogonal bases [21]. Here the 
computational basis a)} is trivial, while equal 

real superpositions can be obtained by manipulation of 
a squeezed vacuum. As shown in Fig. [l] this state is 
produced using an optical parametric amplifier (OPA), 
pumped in the collinear regime. For sufficiently weak in- 
tensities, the OPA generates a good approximation of 
an even superposition |+)=A/+ (|a) + |— a)), where A/+ 
is a normalizing factor p!5l. The odd superposition 
\ — )=JV- {\a) — \—a)) results from the application of the 
gate [17 . In our experiment, this is implemented by the 
beam splitter BSq, and the avalanche photodiode APDq. 
These two states can be used to test the behavior of a 
second gate constituted by BSi and APDi. 

Although the states we can produce do have marked 
signatures of the ideal behavior, such as a change in the 
parity, current technology limits us to some approxima- 
tions of coherent-state qubits [20 , 22 , which are not fully 
suitable for a direct characterization of the gate. In a di- 
rect comparison of the observed output to the ideal, it 
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Figure 1: Left : scheme of the experiment. The main laser is 
a pulsed Ti: Sapphire at A=850 nm, pulsed at a repetition rate 
800 kHz, with a pulse width of 200 fs. The light pulses are 
frequency doubled by a crystal of potassium niobate (SHG) 
in order to produce a pump beam for the optical paramet- 
ric amplifier (OPA). This OPA produces a squeezed vacuum, 
which is then transformed into a squeezed one-photon by pho- 
ton subtraction J^. This is performed by directing a small 
fraction of the beam on the avalanche photodiode APDq by 
a beamsplitter BSo, and conditioning the measurement on a 
detection event. Either the squeezed vacuum or the squeezed 
photon are used as the input of the tested gate, constituted 
itself by a beamsplitter BSi and the APDi . Both beamsplit- 
ters are realized by a sequence of a half- wave plate and a po- 
larizer, so that we could tune both reflectivities to R = 10%. 
This also allow us to take them out of the path of the input 
for a direct measurement. A small portion of the main laser is 
used as a local oscillator (LO) for a homodyne detector (HD); 
the phase cp is scanned by means of a piezo actuator. 
Right : Visual representation of the phenomenological model. 
Symbols are defined in the main text. 



would be hard to deconvolute the errors due to the im- 
perfections of the gate, and the imperfections of the input 
state itself. 

These limitations can be relaxed if one concedes to ex- 
perimentalists some a priori knowledge about the physics 
of their gate : a model can then be derived, and used for 
such a deconvolution, and thus to obtain the behavior 
of the device for ideal inputs. This is somehow similar 
to the analysis in [24 to estimate the origin of decoher- 
ence in polarization quantum gates. In our case, we can 
rely on a simple phenomenological model [22 , which is 
grounded on a more rigorous multi-mode treatment [23 . 

This model is illustrated in Fig. [T] for the case of single 
photon subtraction. We consider a perfect squeezer, i.e. 
a collinear OPA, which can reduce noise by s=e~'^^ times 
the shot noise. The process is spoiled by parasite ampli- 
fication that can be described by a fictive non-collinear 
OPA injected on one side with the squeezed vacuum; its 
gain is expressed as h=cosh.{jr). Once we trace out the 
fictive idler mode, we obtain the following expression for 
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Figure 2: Quadrature distributions of a photon-subtracted 
squeezed photon. The points are the measured values, while 
the solid lines are the prediction of the model, based on the 
outcomes from a squeezed state. Each histogram is calculated 
from 10000 quadrature points. We can attribute some noise 
to fluctuations during the long acquisition time, due to the 
low two-photon subtraction rate (^6 events- s ~"^). 



This is the squeezed state that is fed in the gate. 
This consists of a beam-splitter with transmissivity Tc^l 
(T=0.9 in our implementation), and an avalanche photo- 
diode (APD) of efficiency a^, whose signal is used as a trig- 
ger event. The mode of the squeezed vacuum is filtered 
by a single mode fiber and a spectral filter; this opera- 
tion, though, is performed with finite accuracy. There 
is then a fraction ^ of the events that originates from 
photon subtraction on the correct mode; the remaining 
1—^ are actually non-correlated events that corrupt the 
functioning of the gate. 

The finite efficiency of the homodyne detection (HD) 
has to be taken into account; for the sake of simplicity, we 
introduce the homodyne loss 1—r] before the beamsplitter 
rather than before the HD, then we correct the efficiency 
of the APD as hc/t]. This has the advantage to simplify 
the computations and is strictly equivalent to the use of 
homodyne loss just before the HD and an APD efficiency 
of Hi. In the limit of small efficiency, we can describe the 
action of the APD simply as the destruction operator a 
acting on the reflected mode, that is finally traced out. 

We obtain thus a parametrization of the Wigner func- 
tion, which depends on a few parameters; some of them 
are inherent to the state - the squeezing 5, the parasite 
gain h - other describe the action of the gate - the modal 
purity ^, the transmissivity T - and finally the homodyne 
efficiency r] describes the detection. These can be either 
directly measured by a classical signal, as it is the case for 
T, or inferred by a best fit on the quadrature histograms. 
For our analysis, we extrapolate the results for the gen- 
erated state, by imposing r]^l. The same modeling is 
applied when considering the double photon subtraction. 

For our purpose, it is crucial to note that the function- 
ing of the gate is well described by only two parameters: 
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Figure 3: Fidelities of the output states for arbitrary ideal inputs, a) for our phase gate; b) for a perfect modal purity ^ = 1; c) 
for a device with T=0.99; d) for our device with a=1.2. The angles and are the spherical coordinates on the Bloch sphere. 



the reflectivity of the beamsplitter T, and the modal pu- 
rity ^. These parameters can be accessed experimentally 
by using a set of states as a probes. In our case, coherent 
inputs are not useful, as they are less sensitive to faulty 
events. We thus use the squeezed vacuum, approximating 
the even cat state, as a probe to estimate the parame- 
ters. Then the squeezed photon, approximating the odd 
cat state, is used to check the consistency of the model. 

More in detail, we performed homodyne measurements 
on the squeezed vacuum, from which we estimated its 
density matrix po. This is the starting point of our model, 
from which we can calculate the expected action of the 
gate, given this state as the input. From a fit, we can 
derive the value of ^. At this point, we use a squeezed 
photon as the input; as before its density matrix pi is 
reconstructed by homodyning. The model is now used 
to estimate the output, fixing T and ^ at the same val- 
ues as we found in the first case. Finally, we check the 
consistency of the expected and measured outputs. 

Typical results are shown in Fig.|2] we plot six different 
histograms of quadrature distributions, where the points 
show the experimental results, and the solid lines the 
prediction of our model. We remark that the parameters 
are not fitted on the data, but derived from a previous 
measurements on the squeezed vacuum. This gives us 
evidence of the reliability of our predictions. 

By using this model, we can estimate the action 
of the gate on an ideal arbitrary input in the form 
= A/6>,0 (cos + e*"^ sin ||— a)), with and ^ 
defining the Bloch sphere of the qubit [20 . For the am- 
plitude |a| of the coherent states we need an estimation 
of what experimental limits are. For this, we have taken 
the value |a|=0.92 giving maximal fidelity between our 
squeezed photon and an ideal odd superposition. 

Our analysis is summarized in Fig. [3j where we plot 
a table of the fidelity F^^^ of the output state with the 
ideal The fidelities for our experimental condi- 

tions are depicted in Fig. |3^. The first remark concerns 



the directionality of the device. Better results at the poles 
with respect to the equator are expected, due to the fact 
that the high transmission make coherent states largely 
insensitive to the modal purity ^. Furthermore, our anal- 
ysis shows that odd superpositions work better than the 
even ones. The limit ^=1 (Fig. [sJd) reveals that it is 
an intrinsic feature of the gate: it results from the facts 
that our APD has a low efficiency and to a lesser extent 
that it can't resolve the photon number. However, as the 
transmissivity T increases (Fig. [sj^), Fq^^I) also increases 
and becomes more regular with respect to 0. Indeed, the 
probability to have several photons reflected to the APD 
decreases when T increases, and therefore the fact that 
the APD does not resolve the number of photons becomes 
less important. For superpositions with ^=7r/4 the limit 
value of 0.83 (when T^l) is almost reached for T=0.99. 
The downside is that the probability of success of the 
gate decreases proportionally to (1— T). Finally, we see 
that when a increases (Fig. |3]i), Fg*,^ also becomes more 
regular with respect to ^, but it reaches lower values, as 
the probability to have several photons sent to the APD 
is higher. 

As a more general comment, we conclude that, al- 
though the value ^=0.83 might be considered satisfac- 
tory in our experiment, the operation of the gate is heav- 
ily influenced by the modal purity. This sets a strong 
requirement for realizing such gates. 

The results in Fig. [3] provide us with extensive in- 
formation about our device, but fail in delivering us 
a conclusive answer on how good is the gate overall. 
As mentioned before, associating a quantum process to 
this operation is non trivial, so we can not formally use 
Jamiolkowski's isomorphism and derive a process matrix 
[25] . However, we can still retain the underlying physical 
idea behind the isomorphism. 

Let us briefly recall Jamiolkowsi's construction; given 
a process f on a space of dimension (i, we can asso- 
ciate univocally a matrix by considering a maximal en- 
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tangled state \'^) = ^J2^=i ^ given choice of 

basis {\j)}i<j<d of the system. The process matrix is 
then obtained through the application of the process £ 
to one of the modes of the entangled state, what can be 
formally written as x=^ ^ ^ ^^^^ ^ state in 

a larger Hilbert space. When a distance between two 
processes is needed, we can then rely on the measures 
holding for states in the extended space; in particular, 
the fidelity between two processes is given by the fidelity 
between the corresponding states: F(f i, f2)=^(xi7 X2) 

This is in some sense an application of quantum paral- 
lelism. Since we are interested in the overall behavior of 
the gate, we need to estimate its action on all the inputs 
at the same time: this amounts to feed in the gate half 
of an entangled pair. In light of these considerations, we 
can adopt as a reasonable figure of merit a fidelity be- 
tween entangled states. Here we will use as a probe the 
state |<|)+) = -^ + but actually the result 

does not depend on this choice (see Appendix). 

We then consider the fidelity F between the entangled 
output as would be produced by our device, and the ideal 
target state \Q)=i ^ a\^^) /\\I ^ a\^^)\\. Our results are 
shown in Fig. |4] the limiting factor of the performances 
is mostly the worst-case scenario of an even superposi- 
tion. We find F=0.78 for the experimentally observed 
gate (C=0.83). 

Let us emphasize that the ideal target state differs from 
the bit fiip state |^+) = ^ (|+)|-) + |-)|+)) which would 
be obtained by an ideal phase gate. This is due to the 
fact that the states a) have different normaliza- 

tion coefficients A/± . Such a problem becomes negligible 
for large enough values of a, i.e. when the two states 
|±a) become nearly orthogonal. This requirement can 
be made quantitative by evaluating the fidelity between 
the target state \Q) and |^^), which is simply : 



tanha^) 



(1) 



This fidelity is plotted on Fig. |5) and gives an idea of the 
required values of |a| for the protocol to work correctly, 
in the sense that cat-state qubit is "good enough" for 
correctly implementing the desired quantum phase gate. 
For our value a=0.92, J'=0.967. 

Summarizing, we have quantified the quality of a sim- 
ple (tt phase shift) quantum gate for "cat-state" qubits, 
by evaluating either the fidelity between an arbitrary 
input state |?/^6>,0) = A/6i,0 (cos + e*"^ sin ||— a)) and 
the corresponding output state, or the action of the gate 
on half an entangled state. This second approach gives 
a single figure, F=0.78 in our case. The method is rela- 
tively independent of the quality of the quantum states 
used to probe the gate, but it requires some a-priori 
knowledge and modeling of the gate. It appears thus 
as a stable approach for obtaining both a "quality eval- 
uation" from a single number, and a full description by 
concatenating simple but efficient theoretical modeling of 
the gate. 
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Figure 4: Fidelity of an entangled output from the phase gate 
with the ideal entangled state: we can notice a linear increase 
of the overlap with the modal purity ^. The red dot corre- 
sponds to our experimental parameters. 




Figure 5: Fidelity between |Q) and the ideal bit-flipped state 
1^^). The red dot corresponds to our experimental parame- 
ters. 
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Appendix : Invariance of the fidelity with the input 
entangled state 

In the main text, we have considered the case when 
a hah of the entangled state |^+) = -^ (|+)|+) + |-)|-)) 
is used as the input of the gate. Here we show that 
the result does not actually depend on the choice of the 
Bell state. For this purpose, let us consider the general 

where corresponds 

corresponds to the op- 



case 



l*) = ;)5(l+)lM> + e'1 
either to |+) or | — ), whereas 
posite ket: |i^) = |— /i). 
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The superoperator £ describing the action of the gate 
can be decomposed in two parts S'^sood) g{bad) 
responding respectively of the correct and faulty events: 
S(9ood)_^(^i _ ^) gibad)^ YiTst, we notice that by lin- 
earity we can calculate separately the action of the gate 
for the correct heralding events and for the faulty ones, 
and then sum the results with the correct weighting. 

Let us now focus on Sidood)^ 

reasoning; be- 
ing similar for ^(^^^). For a given initial state 
P=^x,y=+-^xy\x){y\, the operator ^(fi'^^^) is non lin- 
ear. However, as it is modeled with linear operators 
and partial traces, its non linearity comes only from 
the final normalization. One can thus write it as a 
linear operator E^^ood) giy^j^g a non normalized output 
state £^^^^^\p)^ which is then normalized by its trace: 

S{90od)^p^^^g{good)^pyrj.^^g{good)^p^Y dcuotc with 

the symbol Cxy the action of E^^ood) ^^ie operator 
\x){y\. Under the operator I^S^^ood) ^-^^ g^^^^ |^)(^| is 
transformed into: 



x = ^(I+)(+I<^Cmm + I-)(-I<^C.+ 

e'1-)(+l®C.M + e"'1+){-|®CM.) 
Then, we notice that a|/i) = a'^\v). Introducing 



Cfiu = -jfi^ the ideal target state is: 



/(g)a|^) 
1 



(3) 



{V^\+)W)+e"^y^\-)\n)) (4) 



We then consider the fidelity between the normalized 
state x=x/Tr{x} and This has the explicit expres- 
sion: 



(^Ixl^) 



2Tr{x}(c^.+c,^) 



(5) 



r 



which does not depend on the choice of (/) and is invariant 
by exchanging /i and u. The same reasoning applies when 
considering faulty trigger events. In case one would prefer 



to use the bit-fiipped state as the target state, the fidelity 
would be different, but still independent on the choice of 
the initial Bell state. 
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